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Abstract 

In this paper, we propose two subspace-projection-based precoding schemes, namely, full-projection 
(FP)- and partial-projection (PP)-based precoding, for a cognitive radio multiple-input multiple-output 
(CR-MIMO) network to mitigate its interference to a primary time-division-duplexing (TDD) system. The 
proposed precoding schemes are capable of estimating interference channels between CR and primary 
networks, and incorporating the interference from the primary to the CR system into CR precoding via a 
novel sensing approach. Then, the CR performance and resulting interference of the proposed precoding 
schemes are analyzed and evaluated. By fully projecting the CR transmission onto a null space of the 
interference channels, the FP-based precoding scheme can effectively avoid interfering the primary system 
with boosted CR throughput. While, the PP-based scheme is able to further improve the CR throughput 
by partially projecting its transmission onto the null space. 

Index Terms - Cognitive radio, interference mitigation, MIMO, precoding. 
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I. Introduction 

A cognitive radio (CR) Q-H system may coexist with a primary network on an either 
interference-free or interference-tolerant basis [4]. For the former case, the CR system only exploits 
the unused spectra of the primary network. While, for the latter case, the CR system is allowed 
to share the spectra assigned to primary network under the condition that it must not impose 
detrimental interference on the primary network. Therefore, the interference from the CR network 
to the primary system should be carefully managed and cancelled in order to protect the operation 
of the primary system. Various interference mitigation (IM) techniques applicable to CR networks 
have been reported in [IS). As for multiple-antenna CR networks, transmit beamforming (for single- 
data-stream transmission) [l6l- [fT0l or precoding (for multiple-data-stream transmission) [fTT|-[fT3ll 
is an effective approach to proactively cancel interference from CR transmitters to the primary 
network. On one hand, it steers the CR transmission to avoid interfering with the primary network. 
On the other hand, it exploits the diversity or the multiplexing gain of the multiple-antenna CR 
system to enhance the reliability or efficiency of the CR system. 

However, in the works [l6l- [fT3l . perfect or partial CSI of CR interference channels to primary 
network (CR-primary interference channels) is required at the CR transmitter side to guarantee 
no/constrained interference to the primary system. Therefore, extra signaling between primary 
and CR networks is inevitable to obtain the CSI, which jeopardizes the applicability of these 
beamforming and precoding schemes. A more practical precoding scheme - sensing-projection 
(SP)-based precoding, which learns the CSI using subspace estimation and does not require a 
priori CSI, has been proposed for a CR multiple-input multiple-output (MIMO) transmitter-receiver 
pair coexisting with a primary time-division-duplexing (TDD) system in ifTSll . lfT6ll . However, 
such precoding scheme does not account for the interference from primary transmitters to the 
CR receiver (primary-CR interference), which leads to a CR throughput loss. In [fTTII and lITSll . 
it is proposed to remove the primary-CR interference at the CR receiver via null-space receiver 
beamforming, which sacrifices the CR throughput as well. Moreover, the CR network in [fTTl . [fTSll 
has to work in a TDD mode aligned with the primary system in order to facilitate the null-space 
receiver beamforming. 

In this paper, two enhanced SP-based precoding schemes, namely, full-projection (FP)- and 
partial-projection (PP)-based precoding, are proposed for CR MIMO systems by incorporating the 
primary-CR interference. As the name suggests, the FP-based scheme nulls the CR transmission 
by fully projecting the transmission onto the estimated null space of the CR-primary interference 
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channels. Instead of removing the primary-CR interference using null-space receiver beamforming, 
the proposed precoding schemes account for the primary-CR interference via sensing. This, on one 
hand, improves the CR throughput, and on the other hand, introduces more flexibility into the CR 
deployment, i.e., the CR network does not have to work in a TDD mode as in ifTTl . ifTSll . The PP- 
based precoding can further improve the CR throughput by projecting the CR transmission onto 
a subspace that partially spans the estimated null space of the CR-primary interference channels. 
As a result, the CR throughput is further improved at the cost of introducing extra interference to 
the primary network. 

The remainder of this paper is organised as follows. The system model is given in Section II. 
The principle of the SP-based precoding is introduced in Section III. The new precoding schemes 
are proposed in Section IV. The performance of the proposed precoding schemes is evaluated in 
Section V. Finally, we conclude the paper in Section VI. 

II. System Model and Problem Formulation 

We consider a CR system shown in Fig. 1, where a CR transmitter-receiver pair shares the 
same spectrum with a primary TDD network. Multiple antennas are mounted at the CR nodes 
and possibly at each of the primary users. The CR transmitter, CR receiver, primary base station 
(BS) and the kih primary user are equipped with M^, M^, M55 and Mj. (k = 1, • - • , K) antennas, 
respectively. Block-fading channels are assumed for the primary and CR systems. 

For a narrowband transmission, the received symbol at the CR receiver can be expressed as 

y = HFs + n + z (1) 
where y G C^"^'^ is the received signal vector at the CR receiver, s G C^*""^ and F G C^'""^' 

are the transmit information vector with E{ss^} = I and precoding matrix of the CR transmitter, 

respectively, H G C^^^^* is the channel matrix from the CR transmitter to CR receiver, whose 

elements are independent and identically distributed (i.i.d.) complex Gaussian random variables 

with zero mean and variance aj^, and n G C^"^^^ stands for the additive white Gaussian noise 

(AWGN) vector with zero mean and covariance matrix E{nn^} = a^I. Moreover, z G C^^^^ 

denotes the interference from the primary network to CR receiver. It can be expressed as 

{HurX^, during primary uplink 
(2) 
HdrX^, during primary downlink. 

The overall CR system is shown in Fig. 1. The matrices H^r G C^-^'^k^i^k Q and H^t G 
xEf=i (^^QQ Pig represent the interference channels from all the K active primary users to 
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CR receiver and to CR transmitter, respectively, during primary uplink. Similarly, H^r G C^rxM^s 
in (|2]) together with H^t G C^^^^ts (^^^ pjg stand for the interference matrices from the 
primary BS to CR receiver and to CR transmitter during primary downlink. All these interference 
matrices (Hur, H^t, H^r and H^t) have i.i.d. complex Gaussian random elements with zero mean 
and covariances a^^., a^^, a^^ and a^^, respectively. Moreover, G C^^^^^^^^ and G C^^^^-^ 
are the transmitted signal vectors of all the K primary users and primary BS, respectively. We 
define the covariance matrix of the interference in (|2]) as Z = E{zz^}. 

III. Principle of SP-based Precoding 
The precoding problem for CR transmission can be expressed as the following optimisation 



problem UM / hFF^H^ 



max logodet I H (3) 

F V / 

subject to Tr{FF^} < P^r (4) 
Tr{G,FF^Gf } < r, fc = l,...,L. (5) 

In ([5]), Gk G C^kxMt i-j^^ channel matrix from the CR transmitter to the fcth primary user. Thus, 
the channel matrix from the CR transmitter to all primary users becomes H;^ = [Gf , • • • , G^]^, 
where channel reciprocity is assumed. The constraints on the CR transmission power and the max- 
imum allowed interference perceived at each primary user are given by @ and ([5]), respectively. 

The projected channel singular value decomposition (SVD) or P-SVD precoding has been 
proposed in [fT2| as a suboptimal solution for the optimisation problem ©-(IS]). It can be expressed 

as 1 

¥ = \J^[{^A-Kl'yY (6) 

where = max(0, •), ji denotes the power level for a water-filling (WF) algorithm, and 
and originate from the SVD of the effective CR channel matrix 

H^^H(I-UcU^). (7) 

Its SVD is expressed as = V^A^^Uf , and Uc in (0) is from another SVD Hf^ = 
VgA^^^U^, which is estimated via sensing in the SP precoding [fTSll . [fT6^] as shown in Fig. 2. 

By analogy with the multiple signal classification technique [14], the signal covariance matrix 
is decomposed into signal and noise subspaces to estimate U^, which can be mathematically 
expressed as 
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^-^ = ^E^-^(^)^^tW (8) 

= UAU^ (9) 
= UgAgUc + U.A^U^. (10) 
In ([8]), r^t{i) = Hutx^(i) + n{i) is the iih received symbol at the CR transmitter, and R^t denotes 
the average covariance matrix of the received symbols. An eigenvalue decomposition is then 
performed on R^t in ([9]), where A = diag(Ai, • • • , AMt) is a diagonal matrix with descendingly 
ordered eigenvalues of R^t and U G C^*^^* contains the corresponding eigenvectors. The matrix 
Rut is further decomposed into interference and noise components in (fTOl) with and being 
the first Kp = rank(Hut) and the remaining (Mf — Kp) columns of U, respectively, and A^ and 
A^ being their corresponding eigenvalue matrices. 

IV. New Precoding Schemes 

In this section, we elaborate the CR precoding during the primary downlink. A similar precoding 

for the primary uplink can be easily obtained, which is ignored here for brevity. When incorporating 

the primary-CR interference, the precoding problem for the CR transmitter during the primary 

downlink can be expressed as the following optimisation problem 

/ HFF^H^\ 

max logsdet I + —— — -— (11) 



subject to Tr{FF^} < P^r (12) 
Tr{G,FF^Gf } < r,, fc = l,...,L. (13) 

Then, the precoding matrix for CR transmission during the downlink can be written aJl] 

F, = U4(//J-A^i)+]^ (14) 

where /x^ is the power level for the WF algorithm similar to that in (|6]) and Ud is obtained through 
the following eigenvalue decomposition 

U,A,U^ = Hf(Z + a2l)-iH^ 

= (I - UcUg)^H^(Z + a2l)-iH(I - UcUg). (15) 

It can be seen from ([T4l) and ([T5]) that in order to obtain the CR precoding matrix, the interference- 
plus-noise covariance matrix R^r = Z + a^I needs to be estimated at the CR receiver, besides the 
estimation of the interference subspace U^U^ at the CR transmitter. 

^See our conference contribution [|20l. 
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A. Full-projection-based precoding 

To enable the estimation of U^U^ and R^r, we propose an enhanced precoding scheme. The 
system diagram for this scheme is demonstrated in Fig. 2. Each CR cycle consists of sensing and 
transmission phases. We name the CR transmission during the primary downlink transmission as 
Tl and uplink transmission as T2. For Tl, the space U^U^ is estimated at the CR transmitter 
during the primary uplink according to d8T)-(fT0l) over Lsi symbols. The estimation of R^r is 
performed at the CR receiver at the beginning of the primary downlink for a batch of Ls2 symbols 
via a procedure similar to dSJ. After obtaining these two estimates, the CR transmitter starts 
transmission Tl using the precoding matrix obtained by ([T4l) . Then T2 follows immediately after 
Tl but right before the sensing phase for the next CR cycle. The CR precoding matrix for T2 
can be obtained by other two sensing sessions concurrent with the sensing phase for Tl. The 
FP-based channel projection (|7]) is employed in this proposed precoding scheme. Therefore, it is 
termed as FP precoding. 

It can be seen from Fig. 2 that the proposed FP precoding scheme shifts the CR cycle of 
the SP precoding rightwards in time. By doing this, several benefits are obtained. Firstly, intro- 
ducing CR receiver sensing phases during both the primary downlink and uplink improves the 
CR instantaneous throughput by incorporating the interference-plus-noise covariance matrix into 
precoding, and consequently improves the CR throughput. Secondly, shifting the CR cycle diverts 
part of the CR transmission from the primary downlink to the uplink which reduces the time that 
primary receivers expose themselves to interference from the CR transmitter. This is beneficial 
to the primary network, since primary users are usually more susceptible to CR interference than 
the primary BS. 

Theoretically, the proposed FP precoding can completely cancel the CR-primary interference if 
there is no error in the interference space estimation (ITO] . However, the IM ability of the proposed 
FP precoding degrades rapidly when the CR interference-to- noise ratio, INR= cr^t/^n' drops below 
a threshold. This is due to the fact that in ([TOl) some components in the noise subspace may swap 
with those in the interference subspace when the noise amplitude is relatively large compared 
to the interference channel gain (Jut- This phenomenon is known di^di subspace swap^ \2l\ . 

For low INR, the interference subspace has a high probability to swap with the noise subspace. 
When a subspace swap happens, ([T5]) can be rewritten as 



^The lower bound on the probability of the subspace swap has been investigated in ll22l and ll23l . 
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U,A,U^ ^ (I - U„Uf )^H^(Z + aH)-'li{I - U„U^ 

= UcUf H^(Z + a2l)-iHUGUf (16) 



which means that the preceding matrix F^^ and H;^ span the same space. Therefore, when the 
CR INR is low the average interference power received at primary users can be expressed as 

ir = E{Tr{Hf,F,F^H„J} oc (17) 

This suggests that the average interference power at primary users is proportional to the channel 
gain between CR and primary users at low CR INR. 

The average CR-primary interference in the large CR INR regime becomes 

/r = E{Tr{H^,U,(^,I - A,-i)+U^H„J (18) 

= E{Tr{H^t(Ud - U,)(^,I - A-i)+(U, - U,)''H„t} (19) 

^ E{Tr{Hf,(X^H,O^N^Ud(/iJ - A^i)+U/N(Hf,X)tH,J} (20) 

= a2p,,E{Tr{H^,(X^H„0^(H^tX)^H„J} (21) 

= ^Tr{Q J (22) 

where ^ is due to the fact that H^^U^ = 0; (|20]) is obtained using the fact that JJ^ — 
U, ^ -(X^H,t)^N^U, for high INR m with X ^ [x,(l), x,(2), • • • ,x,(L^i)], and N ^ 
[n(l),n(2), • • • ,n(L^i)]; (EB follows from the independence of X^H^t and N and E{N^YN} = 
cr^Tr{ Y}I for any constant matrix Y. Note that Q^^ = EjxuXu^} in (l22l) is the transmit covariance 
matrix for the primary user. An interesting fact can be observed from (1221) that at high CR INR 
the average received interference at primary users does not depend on the interference channel 
H;^. The average interference is proportional to the channel noise and inversely proportional 
to the sensing length Lsi. 

B. Partial-projection-based precoding 

To further improve the throughput of the CR link, the CR transmitter may null its transmission 
to a subspace partially spanning the interference space. Therefore, we introduce another precoding 
scheme, namely, the PP precoding. 

The PP precoding works in a similar manner to the above proposed FP precoding except for the 
selection of the interference space. For the downlink CR precoding, the CR transmitter first obtains 
A and U via eigenvalue decomposition in ([9]) during uplink sensing. Then, a subspace U^U^ 
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partially spanning the interference space is obtained by choosing m eigenvectors corresponding 
to the first m largest eigenvalues of A, where m can be determined according to various criteria. 
One candidate criterion is 

^^^^ < n/, (23) 

with Mmin — inin(Mt, Y^^^i Mk). We call rt/d the trivial over dominant interference ratio (TDIR). 
This selection process chooses m dominant interference subchannels to form an estimate of the 
interference space and ignores the other (M^in — tu) trivial ones. When equal power is assigned 
to each CR antenna, rt/d stands for the maximum ratio of the resulting and nullified interference 
to the primary receiver. Finally, substituting the estimated subspace U^U^ for UgU^ in (|7]), 
the precoding matrix for the downlink CR transmission can be obtained via ([T4]) . However, 
we may fail to find a value of m satisfying ([23]) . In this case, the proposed FP precoding is used. 
The uplink CR precoding can be performed in the similar manner as the downlink counterpart. 

The joint probability density function (PDF) of the ordered eigenvalues A = [Ai, A2, • • • , ^M^i^l 
of R,t, with Ai > A2 > • • • > Am^,. >(jlisUM 



jaIAi, A2, • • • , ^M^inj - -WM~:h \ — p — ' — p — ' • • • ' p I y^^) 

ip \ J^p J^p J^p J 

where Pp is the transmission power of each primary user antenna and /^(Ai, A2, • • • 5 ^M^in) with 
Ai > A2 > • • • > AM^i, is given by 



hi^l^ A2, • • • , -^M^in) 



n£r(M^ax-o!n£r(M^in-0! 



(25) 



with Mmax — max(Mt, Xl^i ^k)- Therefore, the probability for the occurrence of (E) is 

Pm= /a(Ai, A2, • • • , ^M^in) dXi dX2-" d\M^,^ (26) 



where S ^ {(Ai, A2, • • • , Am^,J| ([23]) n Ai > A2 > • • • > \m^,^ > ^^}. 

In other words, for the PP precoding scheme the probabilities of using the 'real' PP (m satisfying 
([23]) exists) and using FP are Pm and (1 — Pm)^ respectively. Therefore, the CR transmitter uses 
(1 — Pm) J2k=i + Pm^ and J2k=i ^k dcgrccs of freedom (DoF) for interference mitigation 
in the PP and FP precoding schemes, respectively. Meanwhile, the DoF for CR transmission for 
the PP and FP precoding are — (1 — p^) Xl^Li ^k — Pm^ and Mt — Mk, respectively. 

This means that compared to the proposed FP precoding the PP precoding scheme transfers 
Pm{J2f=i ^k — ^) DoF from interference mitigation to CR transmission, which leads to a higher 
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throughput for the CR hnk. It can be seen from (|24])-(|26]) that in the large INR regime, is 
fixed for a given noise power and Pp, i.e., the probabihty of the 'real' PP and FP does not 
change with the interference channel H;^. Considering the fact from (1221) that at high INRs the 
average interference power of FP is fixed and the average interference power resulting from 
'real' PP I^^ is proportional to the square of the interference channel gain a^^, the overall average 
interference of the PP precoding = Pml^^ + (1 ~ Pm)Ih^ is linearly proportional to for 
large INRs. 

V. Numerical Results & Discussions 

In this section, the performance of the proposed precoding schemes is evaluated via simulations. 
We consider a scenario where a CR MIMO system coexists with a primary TDD system which has 
one 2-antenna BS and two single-antenna users. Each CR node is equipped with four antennas, 
i.e., Mt = Mr = 4, = 2, K = 2 and Mi = Ms = 1. The primary network works as 
a downlink-broadcast and an uplink multiple-access system. The primary BS uses perfect zero- 
forcing beamforming at both the primary downlink and uplink. The transmission power of the CR 
and primary networks is 1. All the results are obtained by averaging over 2000 simulation runs. 

First, we evaluate the throughput of the CR system with the proposed precoding schemes 
over different values of signal-to-noise ratios, SNR= cr|^/cr^. In Fig. 3, the throughputs (mutual 
information in (fTTT) averaged over a CR cycle) of the two proposed precoding schemes are 
compared with that of the SP precoding of [fTSl . [fT6l and the P-SVD precoding with perfect CSI of 
|fT2ll . The system setup is as follows: Lsi Ls2 Lt2 50, Lti 350, = = 1, Per = 1, 
and rt/d = 0.1. It can be seen that the proposed FP/PP precoding schemes lead to higher CR 
throughput than the SP precoding, and the throughput gain becomes larger as the SNR increases. 

Fig. 4 evaluates the impact of CR INR on the CR throughput and the resulting CR-primary 
interference under different precoding schemes. It has the same setup as that of Fig. 3 with 
cr^ = 10~^. By comparing Fig. 4(a) with Fig. 4(b), it can be seen that the proposed FP/PP 
precoding schemes outperform the SP counterpart at low INRs, since they lead to higher CR 
throughput without introducing extra interference. At high INRs, both the proposed FP and SP 
precoding schemes have fixed interference, and there is a fairly good agreement between the 
derived and simulated interference of the FP precoding. Another phenomenon which can be seen 
from Fig. 4(b) is that the interference of the SP precoding is slightly smaller than that of the 
FP precoding. This is due to the fact that the sensing of the SP precoding is longer than the 
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uplink sensing of the FP precoding. Moreover, at high INRs the interference of the proposed PP 
precoding is hnearly proportional to the CR INR, which supports our analysis in Section IV.B. 

VI. Conclusions 

In this paper, two SP-based precoding schemes, namely, FP and PP precoding, have been 
proposed for CR MIMO systems to mitigate the interference to the primary network and improve 
the CR throughput. These two precoding schemes are capable of estimating the CSI of interference 
channels between primary and CR networks and can account for the interference from the primary 
system via a novel sensing approach. Therefore, no extra signaling is required between primary 
and CR systems, which consequently eases the deployment of CR networks. The performance of 
the proposed precoding schemes has been evaluated and compared to that of the existing precoding 
approaches. It has been demonstrated that the FP precoding can boost the CR throughput and does 
not introduce extra interference to the primary system in the low INR regime. The PP precoding 
can further improve the CR throughput if the primary system can tolerate some extra interference. 
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